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Abstract. We report results of a spectrophotometric 
study of SBS 1150+599A and discuss the nature of this 
object based upon our data. Our study shows that 
SBS1150+599A is most probably a planetary nebula lo- 
cated in the Galactic halo and not a cataclysmic variable 
as originally proposed by the authors of the Second Byu- 
rakan Survey from low resolution spectroscopy. We have 
further elaborated on the properties of SBS 1150-I-599A 
(now becoming PN G135.9-I-55.9) with tools used for plan- 
etary nebula analysis. Our photoionization models show 
that, in order to match the observational constraints, the 
oxygen abundance in the nebula is probably extremely 
low, around 1/500 solar, which is one order of magnitude 
lower than the most oxygen-poor planetary nebulae known 
so far. This finding has strong implications on our under- 
standing of the formation of planetary nebulae and of the 
evolution of the Galactic halo. 

Key words: ISM: planetary nebulae: general - ISM: 
planetary nebulae: individual: SBS 1150-I-599A (PN 
G135. 9-1-55. 9), - Galaxies: halos ~ Stars: binaries: sym- 
biotic - Stars: cataclysmic variables 



1. Introduction 

SBS 1150-I-599A pertains to the list of stell ar-like objects 
from the Second Byu rakan Survey (SBS) ( [Balayan 1997 ; 
Stepanian et al. 1999 ). The Second Byurakan Survey is 
a low resolution objective prism survey with a limit- 
ing magnitude of B~19™5 (Stepanian, 1994; Bicay et al 



2000). The observing technique and the selection crite- 
ria for the SBS objects have been described by M arkar- 
ian & Ste panian ( 1983 ), Markarian et al. ( 1987 ), and 
Stepanian ( 1994a , I994b| ). Worth mentioning is that the 
selection criteria of the SBS include the presence of a 
strong UV continuum, emission lines, and/or peculiar en- 
ergy distribution as inferred from the objective prism 
spectra. These criteria have been successful in selecting 
objects such as UVX galaxies (Markarian galaxies), as 



well as a broad class of AGNs (Carrasco et al. 1998; 
^tepanian etal. 1999| ). As a by product, a large number of 
peculiar stars, WDs, composite and emission objects, hot 
subdwarfs and other types of objects have been detected 



Send offprint requests to: G. Tovmassian 
* Special Astrophys. Obs., 357147 Nizhnij Arkhyz, Russia 
** currently Centre Astronomico Hispano-Aleman Apdo 511, 
E-04080 Almeria, Spain 



(Balayan 1997; Stepanian etal. 1999). Among them sev- 
eral new Cataclysmic Variables(CV) were discovered. 

Here, we report on a spectrophotometric study of 
SBS 1150-f599A in order to clarify its nature. This ob- 
ject is located at very high Galactic latitude: In = 135.96, 
bn = 55.94. It was first thought to be a CV based upon the 
appearance of its spectrum. It shows strong hydrogen and 
helium lines superposed upon a hot continuum, a charac- 
teristic patt ern fo r CVs. Later, in a short communication, 
Garnavich (|l999| ) claimed that SBS 1150-t-599A is a sym- 
biotic star or a young planetary nebula. Our study shows 
that this object is not a CV or symbiotic, but probably a 
peculiar and rather old halo Planetary Nebula (PN). 

In Section 2, we present the data acquisition and reduc- 
tion processes, in Section 3 we derive the observational pa- 
rameters for SBS 1150-I-599A. In Section 4 we discuss the 
various hypothesis about the nature of SBS 1150+599A, 
and in Section 5 we perform a photoionization modeling 
of the object, bringing constraints on the oxygen abun- 
dance and temperature of the exciting star. In Section 6 
we discuss some implications of our work. 
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Table 1. The Log of Spectral Observations 



Date UT 


JD 


Telescope + Equip. 


Wavelength 


Number 


Exp. 


Site 








Range 


of spectra 


sec. 




1998 March 20 


2450892 


2.1m, B&Ch+CCD2K 


3600-6200 


10 


900/1200 


SPM 


1998 June 06 


2450971 


2.5m, IDS+CCD TEK3 


4600-5500 


2 


1200 


La Palma 


2000 March 19 


2451622 


2.1, LFOSC+CCD 


4000-7170 


2 


1800 


Cananea 


2000 April 11 


2451645 


2.1m, B&Ch+CCD TEKl 


3700-7500 


3 


1200 


SPM 


2000 April 12 


2451646 


6.0m, UAGS+PMCCD 


3600-8260 


1 


900 


SAO 



2. Data acquisition and reduction 

SBS 1150-I-599A was observed with the 6 m telescope of 
the Special Astrophysical Observatory of the Russian 
Academy of Sciences (SAO RAN) in 1986. The obser- 
vations were carried out with the SP-124 spectrograph 
equipped with a 1024-channcl photon counting system 
(IPCS) (Drabek etal. |1985| ), installed at the NasmithI fo- 
cus. No flux calibration is available. Based upon that low 
resolution spectrum, SBS 1150+599A was classified as a 



CV ( [Balayan 1997[ [Stepanian etal. 19991 ). 

Recently we acquired a new spectrum of the object 
on the 6 m telescope with the UAGS spectrograph and a 
CCD detector. Although we obtained the necessary spec- 
trophotometric standards, the sky conditions were poor 
for correct absolute flux measurements. The spectrum was 
obtained in the 3600-8260 A range with a spectral resolu- 
tion of 4.65 A/pix which, combined with the deployed slit 
width, gives overall resolution of 10 AFWHM. 

We also obtained a set of spectra of SBS 1150-K599A 
with the 2.1m telescope of the Observatorio Astronomico 
Nacional de San Pedro Martir (OAN SPM) on 20 of March 
1998. The Boiler & Chivens spectrograph was set in the 
4150 to 6750A range with a resolution of 1.28 A/pix cor- 
responding to 4.8 AFWHM. The slit was set to 1.8 arcsec. 
The seeing was poor. It was not possible to obtain stan- 
dard stars for calibration during the same night. We used 
an averaged sensitivity curve for the corresponding grat- 
ing and settings of the spectrograph to calibrate the data 
from this night. 

Additional spectra of SBS 1150-I-599A were obtained 
at the 2.1m telescope of OAN SPM on 11 of April 2000. 
This time we used a 300 1/mm grating with the same 
spectrograph in order to observe the object in the A 3700- 
7500 A wavelength range with 8.0AFWHM resolution. 
During these observations, in order to cover a larger frac- 
tion of the object, we opened the slit up to 3.0 arcsec al- 
though the seeing was good and did not exceed 1.5 arcsec. 
Two spectrophotometric standards were observed along 
with the object and an He-Ar lamp was used as in the 
previous observations for wavelength calibration. 

The object was also observed on one occasion with the 
2.5 m Isaac Newton Telescope (INT) of Observatorio de 
La Palma. The IDS spectrograph was centered around the 
up region (4600-5500 A). A spectral resolution of 1.6 A 



was reached as measured from the FWHM of arc lines. A 
Cu-Ar-fCuNe lamp was used for wavelength calibration. 
No flux calibration is available. 

We also acquired a low resolution spectrum of 
SBS 1150-)-599A at the 2.1m telescope at Cananea, 
Mexico in March of 2000. The LFOSC spectrograph (Zick- 
graf et al. 1997 ) was deployed to cover the A 4000-7100 A 
wavelength range with «13A FWHM resolution. The 
long slit was deployed at the entrance of this multi-object 
spectrograph 3 arcsec wide as projected on the sky. 

For none of the observations was the slit oriented along 
the parallactic angle. Usually it was simply in East- West 
direction. For the SPM observations, due to the faintness 
of the object and absence of appropriate slit-viewing de- 
tector, pointing and subsequent guiding were not perfect 
and might cause inaccuracy of flux determinations. 

Most of the observed data were reduced using IRAF 
packages after applying standard procedures for bias and 
flat field corrections. The observations acquired at SAO 
were reduced using the MIDAS system. The MIDAS cos- 
mic rays cleaning procedure was also applied to the rest 
of the data. 

The Journal of observations is presented in Table 
We performed an extensive search in publicly available 
surveys and databases for any additional measurements at 
the location of SBS 1150-)-599A. IRAS has no detection of 
any IR source within several arcminutes, even among re- 
jected faint sources. No emission was detected at short 
wavelengths. Neither ROSAT RASS nor UV surveys show 
any emission within corresponding error boxes. Among ra- 
dio surveys the closest registered source lies 6 arcminutes 
away from the object. 

3. Spectral analysis and observational parameters 

The spectrum of the object displays a blue continuum with 
strong narrow emission lines of Hi, Hell and marginally 
detected [O ili] . No trace of He i lines or lines from singly 
ionized ions of O, N or S are detected. In Fig. |l| the 
low/medium resolution spectra of the SBS 1150-I-599A are 
presented at different epochs. 

The spectrum in the upper panel is the best flux- 
calibrated spectrum obtained at SPM with the wide slit 
(it is in reasonably good agreement with the latest SAO 
and Cananea spectra). In the inset of the upper panel 
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the sum of 7 spectra (total exposure of 8400 sec) from 
previous SPM observations is presented. It reaches a S/N 
ratio of 20 in the continuum around Ha. The faint hues 
of [O III ] A 5007 A and He ii A 5411 A appear clearly in that 
spectrum. 

The shape of the continuum and of the line profiles do 
not change from epoch to epoch. There is a small discrep- 
ancy of absolute fluxes between spectra, but this could 
easily be a result of inadequate observational conditions 
rather than variability. These leaves open the question of 
the variability of the SBS 1150+599A. 

In the lower panels of Fig. |l| the individual profiles 
of H/3 and He ii lines from the higher spectral resolution 
spectrum obtained at INT are shown. This spectrum, ob- 
tained with 1.5 A resolution, demonstrates that the emis- 
sion lines of the object are relatively narrow. The mea- 
sured values of FWHM of H/3 and Heii are 2.0 A and 
1.85 A respectively. 

Another important spectroscopic feature is that the 
emission lines of SBS 1150-I-599A are displaced towards 
short wavelengths. This blue shift shows some dispersion 
from line to line and also from epoch to epoch. 

Before we consider the measured fluxes, it is neces- 
sary to mention that the emission lines are spatially ex- 
tended, which indicates an extended object. The edge to 
edge width of the line in the spatial direction on the INT 
spectrum is w 14 pix, which corresponds to an object size 
of « 10". A similar value is obtained from the SPM spec- 
tra. 

Table ^ lists the parameters of the emission lines. 
Equivalent widths and relative fluxes were measured in 
all the acquired spectra and we decided to present the 
whole range of obtained values, which gives a good idea 
about the errors of measurement. The spread is up to 35% 
in some line fluxes, however relative numbers are usually 
more consistent among the various spectra. We found sys- 
tematically lower values of the Ha/H/3 ratio in spectra 
obtained at SPM. We suspect that there were problems 
with the linearity of the CCD used in those observations 
and that the very intense Ha line was veiled. We favor val- 
ues closer to the upper limit (2.7) in Table The ratios 
of lines comparable in intensity are in better agreement 
from different observations. In the case of absolute fluxes, 
we used the highest values obtained in the best (but not 
perfect) conditions. The same is true of FWHM and ra- 
dial velocity shifts. We only present the values from the 
higher resolution spectra. Generally the resolution is not 
so relevant in the case of line center measurements. But 
we are doubtful about instrumental line proflles and pre- 
fer not to rely on them. It should be noted, however, that 
in all the low resolution spectra the blueshifts of the lines 
are higher. 



4. Nature of the object 

The optical spectra of SBS 1150-I-599A that we have at 
our disposal do not permit us to classify it unambiguously 
into any known type of object. Below we will discuss the 
most plausible options and will seek to identify the nature 
of the object and model it. 



4-.1. A Cataclysmic Variable ? 

Although SBS 1150-I-599A was originally classifled as a 
CV by the authors of the SBS, a careful examination 
of the spectra rejects such an interpretation completely. 
CVs are close binary systems where stellar components 
are white an red dwarfs, and the radiation is dominated 
by the outcome of accretion processes, usually trough disc 
accretion. CV systems are observed mostly in our vicin- 
ity (100-500 pc). The emission lines originate in accretion 
discs, they are broad due to the dispersion of velocities 
in the disc, and show radial velocity variations on short 
terms, because of orbital motion in the binary (Warner 



19951 ). Meanwhile, the emission lines of SBS 1150-I-599A 



are narrow (< 2 A) and they lack any structure or radial 
velocity changes on a nightly basis. Most importantly the 
lines show persistent blueshift of about 3 A at all epochs. 
A CV with emission lines shifted to that extend should be 
a high inclination system, thus in a time span of several 
hours (5-6 hours are half of the longest period of CVs) 
should demonstrate signiflcant radial velocity variations. 

As mentioned above, in CVs one should expect much 
wider line proflles, due to the presence of an accretion disc 
or accretion flow. Accreting matter is the major source of 
light in CVs and even in magnetic CVs, where no disc is 
formed, emission lines should contain much broader com- 



ponents, variable with orbital phase (Schwope |1997 ). 

Finally, CVs with strong He ii emission as a rule con- 
tain also significant neutral helium and show strong emis- 
sion of He I A 4471, 5876 A, though they are not known to 
have any measurable [O ill] A5007 A emission. 

4-2. A Symbiotic Star ? 

The spectrum of SBS 1150-I-599A looks more like that of a 
symbiotic star (SS) than of a CV. SSs are binary systems 
that contain a giant secondary component (at least for 
80% of them, the rest containing a Mira variable) instead 
of a late type dwarf as in CVs. They can be discovered 
at much larger distances, as far as in the halo or in the 
Magellanic Clouds. Quite a few SSs are known to be halo 
objects. They are also known to contain a variety of neb- 
ular lines ([Oiii], [Neiii], [Nev] and [Fevii]) not seen in 
CVs. 

For years, symbiotic stars were a collection of strange 
objects. However starting from 80's, two c atalog s of SSs 
and s uspected objects were published (Allen 1984 ; Kenyon 
1986 ) consolidating a group of characteristics necessary 
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Table 2. Measurements of spectral lines 



Emission 


Equivalent Width 


T7i„-^^a/- 1 ^— 14\ 

rlUX (^XiU ) 


Relative 


r WrlM 


c 

Line Center 


Line 


xiange 


( —2 —1 \ 

(^erg cm s ) 


Flux Range 




oniit yPs. ) 


He 


-3-5 


0.13 


0.11-0.14 


- 


- 


TT r 

Ho 


-6-12 


0.3 


0.23-0.34 






H7 


-16-18 


0.67 


0.41-0.57 






H/3 


-50-80 


1.19 


1.0 


2.0 


-2.83 


Hq 


-400-600 


3.01 


2.46-2.9 






He 11 A 4686A 


-30-44 


1.09 


0.7-0.92 


1.84 


-2.68 


[0 III] A 5007A 


-2-3 


0.037 


0.03-0.035 






Hen A5411A 


-5-7 


0.064 


0.05-0.065 






He I A5876A 






< 0.01 






[Nil] A6584A 






< 0.005 






[Sii] A6725A 






< 0.005 







flux measurements are based upon the highest values obtained from SPM 2000 observations with a wide slit. 
^ as measured on the higher resolution INT spectrum. 

the shifts as measured on the INT spectrum correspond to a heliocentric velocity of -190±2 km/sec. 
the shifts measured on lower resolution SPM spectra are systematically higher: of order of -260 km/sec 



for an object to be classified as a symbiotic star. More 
recently in a new catalog of symbiotic stars (Belczinski 
et al. 2000) the classification criteria of symbiotics were 
refined, based on a compromise between a large variety of 
features displayed by more than 180 objects classified as 
SSs. 

The three basic characteristics that an object should 
fulfill in order to be classified as a SS are: 

— the presence of absorption lines of a late type giant; 

— the presence of emission lines of H i and He i in combi- 
nation either with lines from ions with ionization po- 
tential >35eV like [Oiii], or with a A-F type contin- 
uum with additional shell absorptions in hydrogen and 
helium lines; 

— the presence of the A 6825 A emission feature, even if 
no features of the cool star are observed. 

SBS 1150-t-599A apart from having high ionization 
emission lines does not s atisfy any other criteria estab- 
lished by Belczinski et al. (200C). Neither our spectra, nor 
available IR surveys provide any evidence for the presence 
of late giant in the immediate neighborhood of the object. 

Another argument against a SS classification is that 
similarly to CVs, SSs exhibit complicated line profiles and 
strong variability. A broad component, in addition to the 
narrow component, should be observed at least in forbid- 
den lines. Otherwise, if higher ionization lines are narrow 
and featureless, one should observe central reversals or 
other structures in the Balmer lines (see for example the 
atlas of high resolution line profiles of SSs of Ivison et al. 
(|199|). Yet, the lines of SBS 1150-h599A do not show any 
deviation from a single Gaussian profile or extended wings 
at the base in any of the observed epochs. 



Taking into account these arguments, we find it very 
unlikely that the object in question is a symbiotic. How- 
ever we cannot exclude this possibility completely. The 
uncertainty with the possible radial velocity variability of 
the lines leaves the question open until new extensive data 
can be collected. 



4-. 3. A Local Group Galaxy ? 

The object being blue shifted, if it were a galaxy it would 
necessarily belong to the Local Group, impl ying t hat its 
distance would be < 1.2 Mpc (van den Bergh 1999 ). From 
the observed H/? fiux, we deduce that the total H/3 lu- 
minosity would be < 2 X 10'^ L©. This is much less than 



the typical H/3 luminosities of Hii galaxies, which range 
between lO^^ - 10^^ erg s'^ (Salzer el al. |l989[ ). During 



the phase where O stars are present, an ionizing star clus- 
ter with a total mass of 10^ Mq (typically the mass of a 
globular cluster, Aguilar et al. 1988D produces about 10^^ 
- 10^^ ionizing photons per second, according to stellar 
population synthesis models (Leitherer & Heckman 1995). 
This implies that, if such a cluster were powering our ob- 
ject, the surrounding gas would absorb only about 1/1000 
of the available ionizing photons. There is no sign of pres- 
ence of large quantities of dust in this object (from the ob- 
served continuum and from the Ha/H/3 ratio) that would 
absorb the radiation, so most of the cluster's ionizing ra- 
diation would leak out. Qualitatively, this would explain 
why the observed HellA4686 A/H/3 ratio is so high (in 
Hii galaxies it is at most equal to a few percent, see e.g. 



Izotov & Thuan 1998) and why no line of Hei is seen. 



However, the presence of an He 11 emission line in an H 11 
galaxy is generally associated with a broad H e 11 fe ature 
produced by Wolf-Rayet stars (Izotov & Thuan 1998 ). One 
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Fig. 1. The spectrum of SBS 1150+599A. In the upper panel the single spectrum obtained with a wider slit is presented. 
In the inset of the upper panel the portion of the summed and averaged spectra is shown to outline the presence of 
faint lines. In the lower panels the profiles of He ii A 4686A and II/3 from the higher resolution spectrum are presented. 



may expect to see a narrow He ii emission in a galaxy if it 
is due to t he ion ization of gas by post-AGB stars (cf. Bi- 
nette et al. 1994 ). But in that case, the expected equivalent 
width of the emission lines would be much smaller than 
observed in our object and the color of the stellar popu- 
lation would be red, not blue. In our spectra, we see no 
signature of the presence of old stellar populations like the 
CN feature (A 4150-4214 A), or G band (A 4284-4318 A), or 
Mgl (A 5156-5196 A) that are found even in Hii galaxies 
(Raimann et al. ^000| ). One should note however, that, if 
this object were an Hll galaxy, its metallicity would be 
extremely low in order to produce such a small [O iii]/H/3 
ratio, and this could explain why the stellar metallic fea- 
tures would not be seen. 



Perhaps the most compelling argument against 
SBS 1150-I-599A being an H ii galaxy is the extremely blue 
nature of the continuum. It is much bluer than the contin- 
uum ob served in blue compact galaxies (see e.g. Izotov & 
Thuan ( |l998|) and references therein) . Fig. || shows the ob- 
served continuum of SBS 1150+599A in the log Flux - log 
wavelength plane, with blackbody functions of different 
temperatures superimposed on it. It clearly indicates that 
the continuum has a color temperature above 50 000 K. 
Even active galactic nuclei do not show such a contin- 
uum m the th e A 4000-7000 A wavelength range (Morris 
& Ward |1988| ). 
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Fig. 2. The low resolution spectrum of SBS1150+599A 
obtained at Cananea observatory is presented in order 
to show the steepness of the energy distribution in the 
continuum. Blackbody curves of various temperatures are 
overplotted for comparison. 



4-.4- ^ Halo Planetary Nebula ? 

Planetary nebulae are the product of the evolution of post- 
AGB stars which, after having expelled their outer layers 
in stellar winds, become sufficiently hot to ionize their 
gaseous envelope. While the spectrum of our object is by 
no means similar to that of a typical planetary nebula, 
since it shows only one forbidden line between 4000 and 
7000 A, we may note that there are examples of planetary 
nebulae in which the lines that are lacking in the spec- 
trum of SBS 1150+599 A ar e ver y weak (AUer & Czyzak 



1983| ; Kingsburgh & Barlow |1994[ ). At first sight, our spec- 



tra can be interpreted as pertaining to a planetary neb- 
ula with a hot exciting star (to produce significant He ii 
emission), that is density bounded (so that no Hei, [Nil] 
and [Sll] lines are seen), and of very low (or perhaps ex- 
tremely high) metallicity, so that the [Oiii]/H/3 ratio is 
very small, or in very high ionization state. For a PN lo- 
cated in the Galactic halo, one can reasonably eliminate 
the high metallicity option. Note that observed heliocen- 
tric velocity of 190-260 km/sec (see Table 2) is indeed 
compatible with SBS 1150-I-599A belonging to the Galac- 



tic halo population (Beers & Sommer-Larsen 1995; Beers 
et al. pOOO|). 



The observed expansion velocity determined from the 
HWHM is 40 km sec'^ (taking a FWHM of 2 A and an in- 
strumental width of 1.5 A). This is rather high for a PN if 



one compares with the values listed in the catalogue of ex- 
pansion velocities of Galactic planetary nebulae by Wein- 
berger (1989) and in the paper by Dopita etal. ( |1985D on 



PNe in the Magellanic Clouds, but there are well known 
PNe which do show similar expansion velocities, for ex- 



ample PNe with [WR] type nuclei (Peha et al. |2001| ). 

In the following section, we further elaborate on the 
PN hypothesis, with the tools developed for the analysis 
of this type of objects. 

5. A model for SBS 1150+599A as a planetary 
nebula. 

5.1. Photoionization computations 

Given the small number of observational constraints, we 
adopt a very simple model of a planetary nebula which 
consists of a homogeneous gaseous sphere surrounding a 
hot star radiating as a blackbody. In order to put some 
limits on the stellar temperature and oxygen abundances, 
we construct series of photoionization models in which we 
specify the stellar temperature, T,, the number of pho- 
tons Q emitted per second by the star at energies above 
13.6 eV, the gas density n, and the filling factor e (both 
assumed constant within the nebula) . Each series consists 
of 12 models with metallicities ranging from 10"'* solar 
to 30 solar (the high abundance end is shown only for 
didactic purposes since, as mentioned before, we do not 
expect a high metallicity object in the halo). Oxygen is 
the reference element (with the solar abundance in units 
of logO /H-l-12 being taken equal to 8.93 from Anders & 
Grevesse (|l98^)). 

The abundances of the other elements relative to oxy- 
gen follow the prescription of McGaugh (1991). Each 
model is constructed with the photoioniz ation code 
PHOTO as described in Stasihska & Leitherer (|l996|) . The 



computations are performed starting from the center and 
are stopped when the equivalent width of the 11/3 line, 
EW(H/3), becomes equal to 130 A. This value takes into 
account the fact that the spectra do not cover the object 
completely (we have checked that changing EW(H/3) by 
50% will not change fundamentally our results). We do 
not have to consider all possible combinations of Q, n and 
e since for a fixed effective temperature and a given condi- 
tion imposed on EW(H/3), the nebular ionization structure 
depends only on Qne^. 

Figure |^ shows the results of our models, computed 
with Q = 3 X 10"'^phs~^, a nebular filling factor e = 0.2 
(values that are typical for planetary nebulae), and two 
values oin: 100 cm~'^ (circles) and lO'^cm^'^ (squares), for 
different values of T*. The first row of panels corresponds 
to T, = 80 000 K, the second to 100 000 K, the third 
to 150 000 K, and the fourth to 200 000 K. From left to 
right, the panels show the behavior of He ii A 4686 A/H/3, 
[Oiii] A5007A/H/3, HeiA5876A/H/3, [Nii] A6584A/H/3, 
and Ha/H/J as a function of log(0 /H) -I- 12. The obser- 
vational constraints are shown by a thick horizontal line 
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Fig. 3. Sequences of photoionization models with varying 0/H. Each row of panels corresponds to a different value of 
T* as indicated on the leftmost panel. Circles: models with n— 100 cm~^; squares: models with n= 10*cm~'^. All of 
the models are constructed using Q = 3 x 10'*'^ phs~^, and e = 0.2, and are density bounded in order to give EW(H/3) 
— 130 A (see text). The observed line ratios are indicated by horizontal lines (dashed lines in the case of upper limits, 
thick continuous lines for measured values, thin lines for the error bars. 



if the relevant line is measured (the thin horizontal lines 
show the error bars) or a dashed line if only an upper limit 
is available. 

All the models show qualitatively a similar behavior. 
For a given value of T, , the effect of increasing the density 
at a given O /H is to increase the ionization parameter, so 
that the ratios He i A 5876 A, [Nil] and [Oiii] decrease. As 
O /H increases, cooling becomes more important and the 
electron temperature drops. This produces an increase in 
the He ii A 4686 A/H/3 and He i A 5876 A/H/3 ratios, due to 
the temperature dependence of the emissivities of these 



lines. On the other hand, [O ill] /H/3 and [Nil] /H/3 first 
increase (abundance effect) then decrease (cooling is be- 
ing gradually shifted from optical to infrared lines). At 
very high metallicities, these ratios rise again, because the 
[O III] and [N ii] lines are now mainly produced by recom- 
binations, due to the very low electron temperature. At 
metallicies less than solar, the values of Ha/H/3 in our 
models decrease with increasing O /H , while the recom- 
bination value (Brocklehurst |1971 ) would increase. This is 
due to the fact that, at high electron temperature, coUi- 
sional excitation contributes to the Ha emission. Clearly, 
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our models show that there is quite a range of input pa- 
rameters that are compatible with our spectroscopic ob- 
servational constraints (if one dismisses the Ha/H/3 ra- 
tio which poses a problem to which we will come back 
below). While models with T, — 80 000 K are not com- 
patible with the lack of detection of He i A 5876 A, models 
with Th. between 100 000 K and 200 000 K are consistent 
with the observed constraints in a range of metallicities 
that depend essentially on the value of the ionization pa- 
rameter (which, in our parametrisation of the problem, is 
only dependent upon the density). For each value of n, the 
value of 0/H which matches the observed [Oiii]/H/3 ratio 
increases with T,, since more and more oxygen ions are 
found in ionization stages higher than 0"*"+. At a given 
T*, the solution for 0/H increases with n for the same 
reason. Clearly, one needs more observational constraints 
to determine 0/H and T,. 

Note that, as seen in Fig. ^ , the observed continuum is 
compatible with any value of T* above 70 000 K and does 
not provide additional constraint. 

5.2. The problem of Ha/Hf3 

It is not satisfactory that the ratio produced by the mod- 
els is much higher than observed. Our computations as- 
sume case B. Case A would better apply in the case of 
our object which is extremely density bounded, but the 
effect on Ha/H/3 would remain virtually unchanged. In 
our models, coUisional excitation of Ha strongly increases 
Ha/H/3 with respect to the recombination value. In order 
to avoid coUisional excitation, one needs to lower the elec- 
tron temperature. This can be achieved by changing the 
geometry so as to increase the fraction of trace neutral hy- 
drogen and, thus, enhance H Ly^ cooling. We have tried 
to obtain such models, but they required very unrealistic 
conditions. 

Since we are not absolutely certain about the calibra- 
tion of the Ha/H/3 ratio in our spectra, we intend to ob- 
tain additional spectroscopy of SBS 1150-I-599A in order 
to check the true value of this ratio. We note, however. 



that in the lists of PNe pu blishe d by AUer & Keyes ( |1987D 
or Kingsburgh & Barlow ( 1994 ) , there are other planetary 
nebulae with Ha/H/3 ratios similar to the one we find here 
for SBS 1150-I-599A. The authors of those papers do not 
comment on their finding, except to say that in those cases 
the reddening to the object was considered zero. 

Interestingly, in some respects, our o bject is simi- 
lar to the one found by Skillman et al. ( 1989| ) in the 
dwarf irregular galaxy Leo A, which shows relatively 
strong He II A 4686 A, no detectable [O ii] A 3727 Aand no 
detectable continuum emission and which they interpret 
as being a PN. Their object shows an Ha/H/3 of 2.3, which 
they tentatively attribute to scattering by dust near the 
PN. 

There are astrophysical objects (old novae, cataclysmic 
variables. Be stars) where the Ha/H/3 is notoriously far 



below the recombination value, reaching values of the or- 
der of one. However, the conditions invoked to explain the 
flat observed Balmer decrements, local thermal equilib- 
rium (Williams 1980D or stimulated emission by a strong 
radiation field (Elitzur et al. 1983| ; Williams & Shipman 
1988) require a high density, which is excluded in the case 
of our object, since it shows the forbidden line of [Oiii]. 

5.3. Distance and structural properties of 
SBS1150+599A 

We now make use of other properties of SBS 1150-1-599 A, 
to further constrain the models. From the observed angu- 
lar radius and H/3 flux (assuming that interstellar extinc- 
tion is negligible, which is likely from the observed value 
of Ha/H/3 and from the high Galactic latitude of the ob- 
ject) we can compute the distance of the object assuming 
a certain mass and filling factor of the gas emitting in H/3. 
Since our object is density bounded, such a procedure (the 
Shklovsky method) provides a reasonable estimate of the 
true distance to the object. We have 



0.14 



-0.2 a-0.6 



M 



neb 



where d is the distance in kpc, Fnp is the H/3 flux in 
ergcm^^s"^, is the angular radius in arcsec and Mneb 
is the nebular mass in solar masses. We adopt an angu- 
lar radius of 5", corresponding to the spatial extension 
of the emission lines in our spectra and a total H/3 flux 
Fh/3= 2.4x10"^'* ergcm~^s~^, assuming that the nebula 
is round. 

Taking Mnob = O.2M0 and e = 0.2 gives d = 20.4 kpc 
which is within acceptable distances for a halo object 
(Zaritsky 1999). Taking Mneb = O.IM© rather than 
0.2Mo gives d = 15.5 kpc, taking e = 1 rather than 0.2 
gives d = 14.8 kpc. 

Adopting as a working hypothesis the values Afncb = 
0.2 M0 and e = 0.2, we obtain a nebular radius Rnob of 
0.5 pc and a gas density of 60cm~'^. The values of Rncb 
and n obtained using various combinations of Mneb and e, 
are reported in Table ^. 

It is easily shown that, for a star of given T*, the ob- 
served EW(H/3) can be obtained by a whole family of mod- 
els sharing the same value of {M^chn/Q) (for a given elec- 
tron temperature). From the models we have computed, 
and which return the value of Mneb that corresponds to 
EW(H/3) = 130 A, we can readily infer the value of Q that 
corresponds to a chosen Mneb- These values, which pro- 
duce acceptable models, are reported in Table ||, for dif- 
ferent choices of T^, . One sees that the values of Q are all 
within a factor 5 of the input value in our models pre- 
sented in Fig. ^ The densities reported in Table ^ for 
SBS1150-h599A are of the order of 50cm-3. Therefore, 
of the series of models presented in Fig. ||, those with 
n^lO** cm~^ are not relevant for our object but those with 
n=100cm~^ have an ionization parameter roughly com- 
patible with the results from Table |[ We can then infer 
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Table 3. Derivation of the structural parameters of 
SBS1150+599A 



These numbers are not to be taken too literally, of 
course. Our model is extremely simplistic. Hydrodynam- 
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an approximate oxygen abundance for SBS 1150+599A by 
finding where the sequences represented by circles meet 
the observed value of [Oiii] A5007A/H/3 in Fig. |. These 
abundances are reported in Table ^ below the correspond- 
ing T,. The accepted values still span a large range, but 
oxygen abundances over one fiftieth of solar are clearly 
eliminated. This makes of SBS 1150+599A by far the most 
oxygen-poor PN known. 

By using a grid of stellar evolution tracks for po st-AGB 
stars (we use the one computed by Blocker ( 1995 ) for hy- 
drogen burning PN nuclei), it is possible to go one step 
further, if one assumes that the transition time between 
the PN ejection and the beginning of the evolution of the 
nucleus off the AGB is zero. We first estimate the expan- 
sion time, from the formula toxp=^neb/i'exp- The derived 
number depends on the assumptions on Mneb and e (see 
Table but is always around 10 000 yr. Then, from Fig. 
^, which has been obtained by interpolating the grid of 
Blocker (1995) on a finer grid in star masses, we derive 
M^[Q) and M:,(r*), the values of the PN nuclei masses 
that correspond to Q and T* respectively. These values are 
reported in Table ||. If the central star evolved exactly as 
predicted by Blocker ( |1995| ) and if SBS 1150-I-599A were 
perfectly represented by our toy model, the true central 
star mass would be obtained for the value of T* which 
makes M^,{Q) and Mh<(T<,) equal. From an inspection of 
Table ||, one finds for the central star of SBS 1150-I-599A 
a temperature of ^ 150 000 K and a mass of ^ 0.58- 
0.59 M0. 



ical computations (Bobrowsky & Zipoy 1989; Mellema 
1994) have warned against a blind use of expansion times 
and we know that PNe are seldom spherical with a uni- 
form density. Besides, it has been shown empirically by 
studying samples of planetary nebulae in various galaxies, 
that the central stars of planetary nebulae evolve differ- 
ently according to the g alactic context (Stasihska et al. 



1998)t Dopita et al. |l99'^ ) 



Nevertheless, the fact that, with our simple approach, 
we are able to give a reasonable picture of SBS 1150-I-599A 
is extremely encouraging. As a consequence, we favor for 
SBS 1150+599A an oxygen abundance of log 0/H-f 12 - 
5.8-6.8. 

5.4. The helium abundance in SBS 1150+599A 

The observed Hen A4686A/H/3 leads to a He++/H ra- 
tio of 0.072 ± 0.010. The observed limit on He i A 5876 A 
leads to He+/H < 0.009. Unfortunately, due to the qual- 
ity of our data, these estimates are not sufficiently ac- 
curate to provide any useful indication of the primordial 
helium abundance. The only thing which can be said at 
this stage is that the helium abundance we are finding is 
consistent with our interpretation that SBS 1150+599A is 
an extremely oxygen poor PN. 

5.5. Need for further observations 

One way to better constrain the oxygen abundance 
(and to test the hypothesis of mixing, see Sect. 6) in 
SBS 1150-I-599A would be through deeper spectroscopy in 
the optical, which should reveal the weak [Ne ill] A 3869 A, 
[Nev]A3426A lines, and possibly also the lines from 
the low excitation e.g., [On] A 3727 A, [N 11] A 6584 A. Our 
models show that even the [O ill] A 4363 A line could be 
measurable. This would considerably improve the qual- 
ity of the oxygen abundance determination. Ultraviolet 
spectra would bear important information. For example, 
they would place limits on the carbon abundance and 
on the ionization structure of the nebula through the 
[O iv] A 1394 A line. These lines are expected to be measur- 
able, despite the low metallicity of the object, due to the 
high electron temperature which boosts the coUisionally 
excited ultraviolet lines with respect to the optical ones. 
A direct observational constraint on the nebular density 
would be very valuable too. The optical doublets which 
are used for that ([Sii] A 6716, 6731 A and [A iv] A 4711, 
4740 A) would be probably too weak to be measured, but 
one may hope to measure the [Ciii]A1907, 1909 A doublet 
in the ultraviolet. 

We show in Fig.|^ the predicted intensities relative 
to H/3 of some important lines, for the same models 
as in Fig.^. These predictions are purely indicative and 
are given only to estimate the feasibility of future ob- 
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servations. For example the intensity of [Nev]A3426A 
strongly depends on the assumed spectral energy dis- 
tribution (here a simple blackbody was used). That of 
[Cm] A 1907, 1909 A is probably larger than shown here 
since C/0 ratios in PNe tend to be larger (Henry et al . 
2000| ) than given by the prescription of McGaugh ( 1991 ) 
used in these exploratory models. 

6. Discussion 



Our main conclusion is that, most probably, 
SBS 1150+599A is a halo planetary nebula. With 
the nomenclature of the Strasb ourg-E SO catalogue of 
planetary nebulae (Acker at el. 1992| ) it becomes PN 
G135.9-I-55.9. Our analysis favors an oxygen abundance 
lower than 1/100 and probably around 1/500 of solar. 
This would make it the most oxygen-poor PN known 
so far. Among the dozen of known Galactic halo PNe 



(Conlon etal. 1993; Napiwotzki etal. 1994; Howard etal 



1997 and references therein; Jacoby et al. 1997), the most 



oxygen-poor is K648 with a logO/H-|-12 estimated at 
7.61. 

How does this fit with our knowledge of the formation 
and evolution of planetary nebulae and on the evolution 
of the Galactic halo? 

It is obviously not surprising to find a low metallicity 
object in the halo, hundreds of stars have been discov- 
ered with metallicities (as referred to the iron abundance) 
[Fe/H]< — 2dex, the lowest ones reaching — 4dex (see 



gen abundances measured in stars are around 1/200 so- 
lar (Boesgaard etal. 1999 ; Israelian etal. 2001 ). It should 
be remembered, though, that oxygen abundance measure- 
ments in halo stars are difficu lt an d unc ertain (see the 
discussion by Mc William 1997 Hill 2001 and the papers 
edited by Barbuy 2001 ) . Our determination of the oxygen 
abundance in SBS 1150+599A is not very accurate either, 
and it would be important to constrain it better with ap- 
propriate observations. Indeed, if the oxygen abundance 
in SBS 1150+599 A reflects the oxygen content of the ma- 
terial out of which the progenitor star was formed, as is 
expected from theoretical models of AGB stars (Forestini 
& Charbonnel 1997), a proper determination of this abun- 
dance would give interesting clues to the understanding of 
the formation of the Galactic halo and of the first stars. 



There is, however, the possibility that the oxygen 
abundance in SBS 1150-I-599A has been modified from the 
pristine oxygen abundance by mixing processes in the pro- 
genitor star. Indeed, differences in abundance ratios from 
star to star in globular clusters giants seem to indicate 
that in those stars deep mixing has brought to the surface 
the products of proton capture in the hydrogen burning 
shell. (Ivans etal. 199£; Weiss et al. 2000 and references 



e.g. references in Beers 1999[ Norris 1999). Oxygen abun- etal 



therein) However, this is not seen in field stars (Gratton 
et al. |2000| ; Charbonnel & Palacios ^001 ) . Usually, abun- 
dance ratios observed in PNe do not indicate that oxygen 
has been depleted in the progenitor by mixing (Richer 
1998). There is however the case of the halo plane- 



dances, however, have been determined in only a fraction 
of these metal-poor stars (Laird 1999). The lowest oxy- 



tary nebula BB-1 which shows a Ne /O ratio much h igher 
than the remaining planetary nebulae (Howard et al. 1997 
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Fig. 5. Predicted values of some line ratios for the same models as shown in Fi^. 



and references therein), and for which destruction of 
in the ON cycle has been suggested. 

If the progenitor of SBS 1150+599 A was indeed one of 
the most metal poor objects known in the galaxy, could it 
have produced a planetary nebula that we see now? With 
such a low metallicity, the star must be older than the 
most metal-poor globular clusters, i.e. about 14Myr (Van- 
denBerg 1999D an d from stellar evolution models (Char- 
bonnel etal. 1999) we derive that the initial stellar mass 
would be 0.8 — 1.0 M0. If we adopt the initial-final mass 



relation of Weideman (1987) or even that of Richer etal. 
( pW7| ), the mass of the planetary nebula nucleus would 
be around 0.55 — 0.555 Mq. According to the post- AGE 



stellar evolution models of Blocker (1995 and references 



therein), the evolution of the star off the AGE would be 
so slow that the ejected material would have been totally 



dispersed before being ionized. There are however several 
possibilities to reconcile the low metallicity of the object 
and the fact that it has produced a planetary nebula vis- 
ible now. First, the star could have been formed from 
infalling material, according to a scenario discussed by 
Carney ( 1999 ). Then its age would be smaller than that 
derived for the bulk of metal poor globular clusters, and 
its mass larger than 1 Mq. Second, there are many uncer- 
tainties in the derivation of the initial-final mass relation 



(Weideman |l987l ; Bragaglia et al. |l995t Richer et al. |1997D 
and it is possible that at low metallicities the final mass 
would be larger, for example due to decreased mass loss. 
One cannot exclude some intrinsic scatter in the initial- 
final mass relation, for example due to close binary stel- 
lar evolution. Third, the stellar post-AGB evolution itself 



can be questioned. For example, Kato & Hachisu (1993) 



12 



G.H. Tovmassian et al. 



argue that computations using OPAL opacities drastically 
shorten the evolutionary time scale. 

For all of these reasons, then, our interpretation that 
SBS1150+599A is an extremely oxygen-poor halo PN 
does make sense. Conversely, a more detailed study of this 
object would provide valuable constraints on the evolu- 
tion intermediate mass stars of low metallicity and on the 
chemical evolution of the galaxy. 

Is SBS 1150-1-599 A unique? If we take a total enclosed 
mass for our Galaxy (within 200 kpc) of 1.4x lO^^ Mp, a 



mass to light ratio of 100 in solar units (Zaritsky 1999|) , 
and if we adopt a PN production rate of 50 x 10^* 
PN/LQ(Ciardullo |l995D , we arrive to a total of 7000 PNe 
in the Galaxy. From the observed metallicity distribution 
of metal poor stars (Norris 1999 ), one does not expect 
to find a large number of such objects in the Galactic 
halo. However, the criteria for discovering PNe in objec- 
tive prism or imaging surveys are based on the strength 
of the [Qui] hne. SBS 1150-h599A was probably first dis- 
carded as a planetary nebula because of its very weak 
[Oiii] line and because of the absence of other emission 
lines usually seen in planetary nebulae. In this respect, 
it might prove rewarding to reexamine the spectra of the 
objects declared as cataclysmic variables or emission line 
stars in objective prism surveys. 
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